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Abstract. We investigate the effects of Pauli blocking on the properties of hydrogen 
at high pressures, where recent experiments have shown a transition from insulating 
behavior to metal-like conductivity. Since the Pauli principle prevents multiple 
occupation of electron states (Pauli blocking), atomic states disintegrate subsequently 
at high densities (Mott effect). We calculate the energy shifts due to Pauli blocking and 
discuss the Mott effect solving an effective Schrodinger equation for strongly correlated 
systems. The ionization equilibrium is treated on the basis of a chemical approach. 
Results for the ionization equilibrium and the pressure in the region 4 000 K < T < 
20 000 K are presented. We show that the transition to a highly conducting state is 
softer than found in earlier work. A first order phase transition is observed at T < 6 450 
K, but a diffuse transition appears still up to 20 000 K. 



1. Introduction 

The physical properties of dense hydrogen are a topic of high interest, in particular the 
transition of hydrogen to a highly conducting phase, which is considered a type of Mott 
transition. Here, we will consider the effects of Pauli's exclusion principle on high-density 
hydrogen plasmas. The previous studies of dense hydrogen included several hypothetical 
assumptions about the character of the high-density phase [H EJ, El HI [5j El [7] . While 
metallization of solid hydrogen near T = K has not been clearly verified so far 
for pressures of up to 300 GPa [8], metal-like features have been observed in shock- 
compression [8j El [TO] • Metal-like conductivities have been observed around 140 GPa 
and 3 000 K [9J. Recent experiments were able to reach that region as well and to 
provide detailed information on the equation of state (EoS) and the conductivity in the 
Mbar-region The transition to metal-like behavior changes drastically our present 
understanding of the behavior of hydrogen at ultra-high pressures. In the present paper, 
we will show that the most important effect leading to the destruction of bound states 
is Pauli blocking. Due to the Pauli exclusion principle, the free electrons in the plasma 
cannot penetrate into the interior of atoms and molecules. At high densities this leads 
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to an enormous pressure acting on the neutrals which will finally be ionized. The 
effective energy levels of hydrogen, which strongly depend on density and weakly on 
temperature, are introduced into the thermodynamic functions. In order to calculate 
the ionization/dissociation equilibrium we minimize the free energy with respect to the 
composition. Recently, we derived an expression for the free energy of dense hydrogen 
Fr [12| [T3l [H] in the framework of the chemical picture and calculated the degree of 
ionization a and the degree of dissociation (3 as well as the isothermal EoS, the hugoniots 
and the isentropes. Pauli blocking was taken into account by the concept of excluded 
volume, which is based on the idea of space occupation by atoms and molecules. We 
will show that a more fundamental approach based on an effective Schrodinger equation 
[T5| [16], leads to important modifications of the earlier results. 

2. Effective Schrodinger equation and bound states of pairs 

We focus on the interaction between atoms and free electrons and present a microscopic 
treatment based on the underlying Pauli exclusion principle. In the following we will 
use Rydberg units with m e = 1/2, U = l,e 2 /47reo = 2, so that the binding energies of 
the isolated hydrogen atom are simply E® = —1/n 2 . Embedding the hydrogen atom 
in a plasma environment, the interactions with the medium are treated by an effective 
wave equation J2J [15], [16] 

p 2 Mp) - E v(q)Mp + ?) + E HPl (p> q)Mp + q) = EuMp) , (i) 

where V(q) = 8n/q 2 denotes the Fourier transform of the Coulomb potential. The 
center of mass motion has been neglected, assuming the adiabatic limit m e /m p <C 1. 
In general, the plasma Hamiltonian H pl (q) will depend also on P and on the energy, if 
dynamical and retardation effects are taken into account. The plasma Hamiltonian will 
shift the energy eigenvalues E n = E® + AE n and will modify the wave functions ip n {p)- 
In particular, due to the plasma interaction the binding energies may merge with the 
continuum so that the bound states disappear, if the influence of the plasma increases 
with increasing density. This dissolution of bound states is called Mott effect and has 
important consequences for the macroscopic properties of the plasma. Let us evaluate 
now the mean-field energy shift of bound states writing the effective Hamiltonian of 
pairs in a plasma as 

H p \p, q) = £ f e (p + q') [V(q)5(q') - V(q')6(q)), (2) 

f e (p) = l/(exp[/5(p 2 /2m — /i)] + 1) is the Fermi distribution. Within first order, the shift 
of the energy eigenvalues is obtained with the unperturbed wave functions <fi n (p) as 

En-E Q n = Y.<Pn(P)V(p'-p)fe(p)MPl-Y,€XP)V(p'-p)fe(plMp)-W 

p,p' p,p' 

The first term in Eq. ([3]) is the Pauli shift which is due to Pauli blocking, and can be 
rewritten by inserting the Schrodinger equation as 

A £Pauli = £ * (p)(p 2 _ E° n )f e (p)Mp)- (4) 
P 
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A simple expression is found in the low-temperature and low-density limit, where 
the Fermi distribution with the normalization Y^ P fe{p) — n e /2 is concentrated near 
p = 0. In the zero temperature limit, we have a Fermi sphere with Fermi momentum 
Pf = (37r 2 n e ) 1 / 3 . The energy shift of the ground state is then 

A £Pauli = 1^ | 0i(o) |2 = 327rne _ (5) 

At intermediate temperatures we may approximate the Fermi distribution by a 
Boltzmann distribution and find 

A£f auli (T) = 32nn e G(T/T ) ~ 32nn e /[l + (77/16) T/T ] , (6) 

where G(x) = {\/x(l + 1/x) — a/7t(1 — x — x 2 /4)[l — erf(l/y / x)] exp(l/x)}/x 7/ ' 2 , see also 
[T7] . and T = 1 Ryd//cs = 157886 K is the ionization temperature. Similarly, the Fock 
term as the second term in Eq. (j3J) can be evaluated as 

roc 1 

AEf ock = -128n e / dp - = -207rn e . (7) 

It compensates partially the Pauli shift. Since the temperature dependence of the Fock 
shift is similar to that of the Pauli shift we may use in first approximation the same 
temperature function as in (jSj) and the total shift is approximated by 

A £Fock + A £Pauli = 12vme G(T/T ). (8) 

The shift is shown in Fig. (TJ dashed line, indicating a rather steep increase of the 
bound state energy with density. Due to phase space occupation, the ground state 
energy is shifted and may merge with the continuum of scattering states, indicating 
the dissolution of bound states. Considering in Eq. (J2]) the continuum part of the 
spectrum describing scattering states, only the Fock shift contributes to the energy 
shift. The lowest energy in the continuum occurs at p = and is shifted by 
AE Fock (p = 0) = - Eg V(q)f e (q) = -4p F /n = -4(3n e /vr) 1 / 3 . However, the two-particle 
continuum state can only be created at the Fermi momentum since all states below that 
are occupied. Thus the continuum of scattering states begins at pf where we have in 
the zero temperature limit the Fock shift 

AE Fock (p F ) = - J2 V{p- p F )f e {p) = -2p F /n = -2(3/ 7 r) 1 /3 n i/3 ) (9) 
p 

shown also by the dashed line in Fig. [TJ Extrapolating the low-density results to higher 
densities, the ground state disappears in first approximation at a density n e ~ 0.015. 
This corresponds to an average distance of ro — 2a b and is below the Mott criterion. 
The Mott condition r ~ as with (47r/3)n e rQ = 1 expresses the idea that atoms are 
destroyed if the mean distance of the electrons is equal or smaller than the Bohr radius. 
There exist many alternative estimates of the binding energy shift. For example, the 
confined atom model [21 US] assumes that the atom is embedded into a sphere with radius 
r . In first approximation this theory gives the shift (in Rydberg units) AEf a = ir 2 r~Q 2 . 
Correspondingly, the binding energy would disappear at r ~ 3a#, i.e., already at a 
much smaller density (see also Fig. [TJ). Better estimates based on numerical solutions 
of the Schrodinger equation give a value of about r ~ 2clb [19]. Our estimate is in the 
same region. 




Figure 1. Density dependence of the effective ground state energy of hydrogen. 
The low-temperature estimate (dash-dotted line) according to Eq. ((8]) is compared 
with the confined atom approach (dotted line). The solid line corresponds to the 
variational approach. We have shown also the lowering of the continuum edge (dashed 
line) according to Eq. ([9|) . 



3. Evaluation of the mean-field energy shift by variational approach 

According to our estimate, the effective binding energy would disappear at n e ~ 10~ 2 . 
This seems to be too early; the reason is that perturbation approximations tend to 
overestimate effects. Better results may be obtained by variational approximations 
to the solution of the effective Schrodinger equation. To apply the Ritz variational 
approach, we have to symmetrize the Hamiltonian in the effective wave equation 
introducing the function ^ n (p) = V>n(p)[l — fe{p)Y 1 ^ 2 - We will consider the zero 
temperature case and use the ansatz corresponding to a variable Bohr radius, 

V ( p] a ) = Sn^a-^il + p 2 /a 2 y 2 . (10) 

Here a is a parameter which characterizes the occupation in the momentum space. 
In a more refined approach we may take into account that no states below the Fermi 
momentum are available. The shift of the binding energy as function of the density is 
shown in Fig. [H Within a better approximation we take the Fermi function in the zero 
temperature limit, f e (p) = ©(pf —p), and evaluate the Pauli blocking shift integrating 
over the wave function ^oip] «)■ We obtain 

4-7T l"PF 

A£f auli « j-^ I dpp 2 [p 2 /a 2 + «(p) 
\ppjo? - l)p F /a 



(2k) 
4a 



71 



+ arctanQop/a) 



(11) 



(l+pl/a 2 ) 2 

We calculate the temperature dependence of the Pauli blocking term by replacing 
the zero-temperature Fermi function in the interaction term by the finite temperature 
distribution. It can be seen that the temperature dependence of the Pauli blocking term 
becomes weak. Even within the variational approach, the densities where the energy 
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levels disappear and consequently full ionization occurs, are evidently still too low to 
explain the observed effects. We include now the Fock term which is of the same order 
as the Pauli blocking term. Even if the Fock term is not of primary importance for 
the disappearance of the bound state energy, it has to be included in the total shift of 
bound and scattering states to be consistent (so-called conserving approximations). In 
the zero temperature limit we get after some transformations the integral 

AE Fock - -- r Pdp [ PF kdkln {p + k) - -A„3 5 + 3p| 

which reproduces in the low-density limit the value (— 207m e ) given above in Eq. (J7J). 
We may estimate the temperature dependence as above in perturbation theory. For 
convenience of the numerical procedure in the later variational calculations of the free 
energy we constructed an interpolation formula between the Boltzmann and the zero 
temperature limits, i.e. Eqs. (6) and (11), taking into account a few points which we 
have evaluated numerically, 

A ^ ai1 = l { tl ritSr+k) + ""-M • (i3) 

This is nearly identical to the asymptotic representation except that we had to introduce 
a fit function c(T) = (G(T) — l)/3 in order to provide the correct derivative at small 
densities. A similar interpolation can be given for the Fock term 

with the fit parameters g = 261.65, k = 60 000, I = 334 369. A comparison of the density 
dependence according to the interpolations introduced above with numerical estimates 
of the integrals is shown in Fig. [2] for T = 5 000 K. The agreement with the data is 
reasonable for this temperature. In the region of interest 5 000 K < T < 15 000 K, the 
temperature dependence is quite weak. The remaining shifts are smaller and will be 
neglected here, see also [TBI fT8l l20[ I2T] . 



4. Ionization equilibrium and thermodynamics in the chemical picture 

We construct the thermodynamic functions of hydrogen by using a chemical approach 
to the free energy which recently was applied to temperatures between 2 000 and 10 000 
K [121 [221 [231 [2H [131 [261 127]. The effects of pressure dissociation, H 2 # 2H, and 
ionization, H e + p, are taken into account so that the transition from a molecular 
fluid at low temperatures and pressures through a partly dissociated, warm fluid at 
medium temperatures of some thousand Kelvin to a fully ionized, hot plasma above 
10 000 K can be explained. 

We will not go into the details of the free energy expression used earlier. Here we 
take the expression Fh as applied so far [T2l [13j [14] and add the contribution due to the 
energy shifts 

F(V, T, N) = V(n a + 2n m )AE(n' , T) + F H , (15) 
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Figure 2. Results of a numerical evaluation of the Pauli shift (dots) and the Fock 
shift (diamonds) for the temperature T = 5 000 K in comparison to the interpolation 
formulae (113|) and (fT4|) (bold lines) and the Boltzmann approximation (thin lines). 



where n a and the number densities of the free atoms and molecules. The shift 

of the atomic levels is approximated by the sum of Pauli and Fock terms. Furthermore, 
we assumed that molecules are simply composites of two atoms so that the shifts are 
additive. We took into account ionisation and dissociation processes. The degrees of 
ionisation and dissociation defined by [12] 

rii „ n a „ 2n r . 



a 



Pa 



ft,, 



(16) 



ni + n a + 2n m ni + n a + 2n m m + n a + 2n m 

are the variational parameters of our problem. The free energy has to be minimized with 
respect to them. We note that ft a is the relative amount of protons bound in atoms and 
j3 m the relative amount bound in molecules. Due to the balance relation for the total 
proton density we have the simplex relation a + ft a + (3 m = 1. Atoms appear only in a 
rather narrow region of the temperature - density plane. The density dependence of the 
degrees is represented in Fig. [31 We calculated the chemical equilibrium by means of 
a numerical variational procedure based on direct minimization of the free energy [25J. 
We prefer here the minimization of the free energy as compared to the Saha approach 
because it finds all existing minima, including those at the boundaries. The transition 
density is between 10 23 and 10 24 protons/cm 3 . This means we are in the region around 
0.5 g/cm 3 for hydrogen and 1 g/cm 3 for deuterium. The transition to full ionization is 
rather soft. In principle, all thermodynamic functions may be calculated from the free 
energy (1151) by derivatives. Examples of our results for the pressure and the relative 
pressure (in relation to the reference pressure of a fully ionized plasma) are shown in 
Fig. HI We see a phase transition of first order below 6 450 K and a diffusive transition, 
defined by wiggles in the relative pressure p/(n p kBT) up to 20 000 K. 
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Figure 4. Pressure isotherms (left) and isotherms of the relative pressure (right). 



5. Discussion and Conclusions 

In order to understand the transition to metal-like conductivity in hydrogen at high 
pressures, we studied the role of Pauli blocking effects. We calculated the Pauli and 
Fock energy shifts by solving effective Schrodinger equations for strongly correlated 
systems. The ionization and dissociation equilibria were treated within a chemical 
approach. We have shown that Pauli blocking effects have a strong influence on the 
ionization equilibria and the character of the transition in the high pressure region. We 
presented explicit calculations of the ionization and dissociation equilibria from low to 
high densities in the region 4000 K < T < 20 000 K. An even more detailed description 
should take into account higher order contributions to the shifts [20J. 

The transitions to highly conducting states occur at densities around 3 ■ fO 23 
protons/cm 3 . The corresponding pressures are in the region of (0.8 — 1.2) • 10 11 Pa, 
i.e. around a Mbar. The first order phase transitions are softer than observed in earlier 
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work. In the new theory, the first order transition appears only at T < 6450 K, however 
a diffuse phase transition detected by wiggles of the relative pressure remains up to 
20 000 K. 

We acknowledge helpful discussions with Wolf D. Kraeft and financial support from 
the DFG funded SFB 652. 
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